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Abstract 

V.I. Arnold has experimentally established that the limit of the statistics of incomplete quotients 
■ of partial continued fractions of quadratic irrationalities coincides with the Gauss-Kuz'min statistics. 

Below we briefly prove this fact for roots of the equation rx 2 + px = q with fixed p and r (r > 0), and 
with random q, q < R, R — » oo. In Section 3 we estimate the sum of incomplete quotients of the period. 
' According to the obtained bound, prior to the passage to the limit, incomplete quotients in average 

are logarithmically small. We also upper estimate the proportion of the "red" numbers among those 
representable as a sum of two squares. 
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OO ; 1 The statement of the main results and their consideration 

The papers, monographs, and reports of V.I. Arnold dedicated to the statistics of periods of continued 



fractions and their multidimensional generalizations (see [T]-[7]) contain a vast number of experimental facts 
and outline prospects for further investigations. The first step in this direction is the proof of the following 
assertion (established by V.I. Arnold): the statistics of incomplete quotients of the continued fraction of 
solutions to a random quadratic equation with integer coefficients, proceeding to the "thermodynamic" limit, 
turns into the Gauss-Kuz'min statistics. Unfortunately, due to circumstances the proof of this assertion 
performed by V.A. Bykovskii and his followers was not published. Below we adduce an elementary proof of 
the one-parameter version of this assertion. 



$_i ' Any number x S K is representable as a continuous fraction in the form 

c3 ; 

x = a H = [ao; Oi, a 2 , . . .], 

ai H 

a 2 + . . . 

where ao 6 Z, a; € N for all i > 1. 

Following Khinchin [5], we denote by E(^'^ ""'^J ) the set of real numbers which satisfy the conditions 

a Sl = At, a S2 — A 2 , . . . , a Sk = A k : 

here, certainly, all Ai and Sj are positive integers, and all Si are different. This set is the union of a countable 
number of intervals [8]. 

Let ^''"'^(i) stand for the indicator function of the corresponding set E: 



A U ...,A 



\ 0, otherwise, 
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Let 

"U:::::i)-jt^ w * (1) 

be the Lebesgue measure of the set of points of the interval [0, 1) which belong to ■EQ^'"''x f! )- 

Let x+(q) "+" stand for a root of the quadratic equation rx 2 + px = q, i.e., x+(q) = — ^r-^, where 
A = p 2 + Arq. In what follows we understand the random choice as the equiprobable sampling from a finite 
set. 

Theorem 1. Let r £ N, p £ Z, q € Ziere r and p are fixed, q is randomly chosen. Let 

Pr P \Ai' '"'a ) s ^ an d for the probability that x + {q) E E\X\'"''a)' Then 

H — ^oo 



( si,- 


• • , Sfe \ 


/ si,- 


• • , Sfe \ 











Remark 1. One can replace the condition q € {1, . . . , R} with that q £ Z, q < R, A > 0; evidently, this 
does not affect the limit probability. 

Using the Kuz'min theorem on the exponential convergence to ln(l + A ^2+2) )/ one can eas ^y obtain 
the following assertion. 

Corollary of Theorem 1. Let r, p, and q be chosen in the same way as in Theorem 1; let the number 
s, s G {l,...,n}, be chosen randomly. Then as R — > oo i/ie probability that a s — A tends to ln(l + 
^ T5J )/ln2 + 0(l/n). 

Remark 2. T/ie probability Y]^—i Pr P {a) l n (considered in the corollary) with fixed R and n — > oo tends 
to the probability that a number chosen randomly from the period of the continued fraction x + (q) equals A. 

The proof of Theorem 1 is based on the following fundamental idea: one can treat integral {T]) not only 
as the Lebesgue integral, but also as the Ricmann integral (considered in the school course on mathematics) . 
A certain special choice of integral sums actually leads to simplified versions of Theorem 1. In a general 
case, it is convenient to perform the proof on the base of the theory of divergent series. 

The idea to apply the integral Riemann sums is not new; it was used implicitly in the proof of a similar 
correlation in the case of rational numbers with a fixed denominator [9]. However, usually one uses another 
technique for this purpose. The asymptotic fairness of the Gauss-Kuz'min statistics in the case of a fixed 
denominator follows from results obtained by L. Heilbronn and J.W. Porter (see [10]). Not long time ago, 
V.A. Bykovskii, M.O. Avdeeva, A.V. Ustinova [11], [12], [13] obtained the limit statistics for finite fractions, 
whose numerators and denominators belong to a sector or to an arbitrary expanding domain. 

Let us explain the (mentioned in Remark 2) connection with the statistics of the period of a continued 
fraction. In accordance with the Lagrange theorem, a continued fraction for the quadratic irrationality (and 
only for it) is periodic. Therefore the continued fraction for x+(g) takes the form 



—p + \Jp 2 + Arq 

— [ a 0i ■ ■ • , a m, a m+l: ■ ■ ■ , Om+T; Qm+lj ■ ■ ■ i "m+Tj ■ 



2r 



]], (2) 



where T = T(r,p,q) is the length of the period of the continued fraction, m — m(r,p,q) is the length of 
the preperiod. Evidently, one can represent the probability that with random q, q < R, a number randomly 
chosen from the period a m +i, . . . , a m +T coincides with A as the limit lim ^"=1 •^ > r P (a) l n ' Unfortunately, 
we did not succeed in proving that this probability tends to ln(l + A ^2+2) )/ m ^ as i? — > oo. 

In papers [2], [3] V.I. Arnold experimentally studies T(l,p, q). He establish that in average this period 
increases proportionally to the square root of the discriminant: 

f - constVA. (3) 
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Let Tg(q) = T(l,Q,q) be the period of the continued fraction of the square root of q. Let T'(Q) = 
S«=i ^o(?)/<3- This kind of an average was studied experimentally by experts in the theory of numbers, 
who stated the conjecture 

f (Q) ~ const VQlog"(Q), where a <0 (4) 

(see [14] and references therein). 

One can easily obtain the following upper bound for T (Q) [15] : 

f'(Q) < const y/Q. 

In [16] (see also [14]) one proves that the left-hand side of this inequality is asymptotically small in comparison 
with the right-hand one, moreover, if the extended Riemann hypothesis is fulfilled, then the order of their 
difference is not less than log((2) log2_e . 

The well-known lower estimate for Tb(Q) significantly differs from the right-hand side of (U). In [TJ] one 
only proves that Tq(Q) > const log(Q). Note that obtaining the bound Tq{Q) > consty / Qlog Q (Q) would 
have allowed one to solve the Gauss question [17] about the slow growth of the number of classes of real 
quadratic fields. 

There is a known bound for the maximal value of To(q). Put 

Lv^J 

D{Q)=Y,t{Q-u 2 ), (5) 

u=l 

where |_-J is the integer part, r(n) is the quantity of divisors of the number n. Formally speaking, the Dirichlct 
theorem (the equality J2u=i T ( u )/Q = InQ + (27 - 1) + 0{1/^Q)) does not imply that D(Q) ~ s/QbxQ, 
because the values of the function r are not uniform. Reasoning more accurately (see [15j and references 
therein), we obtain the bound 

D(Q) = 0(ln 3 (Q)y/Q). (6) 

The mentioned bound for the maximal value of To(q) (recently it was described by V.I. Arnold in [6]) 
has the form 

T (Q) < D(Q). 

This inequality was first obtained by Hickerson in 1973 [18] (see also [E]). Later Western researchers [20] . 
[2T] succeeded in proving the correlation 

T Q (Q)=0(^Q\nQ). (7) 

Earlier mathematicians of the Leningrad scientific school treated it as a direct corollary of the Dirichlet 
formula and results of [19] (see [22]). Moreover, if the extended Riemann hypothesis is fulfilled, then one 
can analogously obtain the double logarithmic bound 

T(r,p, q) = O(VAlnlnA). 

One can construct a continued fraction for a quadratic irrationality, using the algorithm of successive 
approximations by mediants (the Farey fractions). Based on the connection of this algorithm with explicit 
representations of values of an integer- valued quadratic form (see |23j). in Section 3 we obtain the following 
bound for the sum of partial quotients of the quadratic irrationality within one its period. 

Theorem 2. Let rx 2 +px = q be an arbitrary quadratic equation with integer coefficients and real irrational 
roots. In denotations of (0), (0) we have the inequality 

T(r,p,q) 

^2 a m(r, P ,q)+i < /(A/4), where (8) 
i=i 

f(Q) = 2D(Q) + t(Q), if Q is integer; (9) 
/(Q) = 2 r(Q~i 2 / 4), otherwise. (10) 

ie{l,3,...} 
i 2 <4Q 
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In the case of an odd period T(r,p,q) one can improve this bound, namely, divide the right-hand side of 
inequality by two. 

According to the results of computer experiments, for almost all roots of prime numbers in the form 
4fc + 3, as well as the roots of the doubled numbers in the same form, the left-hand side of inequality ([8]) is 
half as large as the right-hand one. 

Problem 5 in the list of problems stated by Arnold in [7J implies the estimation of rate of growth of the 
average value of elements of the period of a continued fraction. To put it more precisely, the problem is to solve 
the alternative between the power rate of growth and the logarithmic one. Comparing inequality (jSJ) with 
formula ([6]) and the results of the numerical tests (|3I4|) . we conclude that the mean value a is logarithmically 
small as against A. 

Note that the geometric method used for the proof of Theorem 2 enables one to easily obtain assertions 
on the palindromicity of continued fractions in cases when p = or r = 1 (see remarks 5 and 6). We are 
going to dedicate a separate paper to the study of problem 12 in [7J. 

Studying periods lengths To(Q), V.I. Arnold stated the problem of the constructive definition of the 
set K of the so-called "red" numbers: 

K = {Q : Q e N,T (Q) is odd}. 

The set if is a subset of the totality M of positive integers representable as a sum of two squares. The 
explicit representation of the set M is well known (see, for example, [IS]). In particular, prime numbers Q 
of the set M obey the formula: Q mod 4 ^ 3. 

Let K n (or M n ) stand for the part of the set K (the set M) which consists of numbers not greater than n. 
According to the Arnold conjecture, the following limit exists: 

,. K n 

lim = c. 

The sets K and M are closely connected with the generalized Pell equation 

x 2 -Qy 2 = -l. (11) 

Theorem 3. 

1. Q G K equation ill]) is solvable in integer numbers. 

2. Q G M O equation ill\) is solvable in rational numbers. 

3. Let P be the set of prime numbers. With sufficiently large n, 

\K n \<\M n \ n 

p-.peP, ^ 

p mod 4 ^1 



4. QeKDP QeP, Qmod4^3 ^ QeAlDP. 

One can easily calculate that 



n a 

■ mod 4 5^1 



1) = 0.64208. 

pi 



This gives a rough bound for c; computer tests performed for all n from 10 5 to 10 6 show that 0.47 < jf- < 
0.48. 

Theorem 3.1 is well known for experts in the theory of numbers (see, for example, |17l Theorem 12]), as 
well as the absence of simple solvability criteria for equation (TTTj) (i. e., for the membership of K). 

The assertion about the sufficiency in Theorem 3.2 is evident. The proof of the necessity is reduced to 
the proof of the following property: if a positive integer is representable as a sum of two squares of rational 
numbers, then it is representable as a sum of two squares of integers. The proof for the case of two squares 
verbatim coincides with that adduced in [24] Appendix to Chapter 4] for the case of three squares. 

Theorem 3.3 follows from two evident corollaries of the previous items of Theorem 3. Apparently, 
Theorem 3.2 implies that the set M contains sets of all numbers in the form AK, 9K, A9K, etc (here we 
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restrict ourselves only with the multiplication by squares of prime numbers, for which the residue of division 
by 4 differs from 1). On the other hand, Theorem 3.1 implies that K does not intersect with these sets. 
To put it more precisely, if Q € K , then Q is neither divisible by 4 nor by all prime p in the form 4fc + 3. 
Otherwise, considering equality (jlip modulo 4 or modulo p, we get the equation x 2 = — 1; it is well known 
that the latter equation is neither solvable modulo 4 nor modulo p. 

Now Theorem 3.4 is reduced to the assertion that all primes in the form 4k + 1 belong to K . This fact 
was proved by Legendre (see |25j . See |17| ) for the proof. He used it in order to solve the inverse problem 
on the representation of such a number as a sum of two squares. 

Note that Theorem 3 is absent in author's article with the same title submitted to the journal "Functional 
Analysis and Its Applications" 

This paper is written under the bright impression of (video) lectures delivered by V.I. Arnold [4], [5], [6]. 
It resulted from the addition to these lectures which, in turn, appeared as a result of the performed numerical 
experiments and the study of remarkable books [9] , [23] , [UJ . 



2 How to apply the condition "accurate to a zero-measure set" in 
a discrete case (the proof of Theorem 1) 

Let us first prove Theorem 1 in the case of the simplest quadratic irrationality ^Jq. 

Lemma 1. Let q be arbitrarily chosen from the set {n 2 + 1, n 2 + 2, . . . , n 2 + 2n}, n € N. Denote by 
Kilu.'.'.'X) the probability that yfq e E^' - '^J. Then 



( si,. 


• • , Sfc \ 


( si,. 


• • j Sk \ 


Ui>- 


..,A k ) 


=M U. 


■.,A k ) 



Proof of Lemma 1. Let us first note that under the assumptions of the lemma [^J = n. Therefore, by 
definition we have 



2n 

Si, • • • ,Sfc 



; i^EtK^-")- (12) 



Ai,...,A k J 2n 

The complement of the set E^''"'"? ) is representable in the form (J E(b' '"'"b )> *' e '' as 

(B l ,...,B k ):B ] ^A ] 

the union of a countable number of intervals. Consequently, the function I^'"'*X 1S countably continuous, 
therefore, by the Lebesgue theorem one can treat integral |T]) not only as the Lebesgue integral, but also 
as the Riemann integral. In order to calculate it, we divide the segment [0, 1) onto 2n equal parts. We 
choose \Jn 2 + i — n as the point, where we calculate the integrand on the interval [^-, j~) (one can easily 
make sure that < \J n 2 + i — n < 4- for all i = 1, ... , 2n). Then the value of the Riemann integral sum 
coincides with the right-hand side of equality (jT3J) . Lemma 1 is proved. 

Now the proof of Theorem 1 in the case r = 1, p = follows from several remarks. 

Remark 3. For an arbitrary sequence s n and any sequence of nonnegative numbers t n such that Y^n=i = 
oo, the convergence s n — > [i implies that 

tlS! + . . . + t n S n 

Sn = — — — > H. (13) 

ti + . . . + t n 

In the theory of divergent series this property is known as the regularity of the Riesz summation method \2b\ 
Chap. 3, Theorem 12]. 

Remark 4. The assertion of Theorem 1 on the convergence of a sequence P n to follows from the conver- 
gence of a subsequence P n to a, if lim — ► 1. 

i— >oo 

Really, it is evident that 

n i n i+l 

~^Pn z < Pr < —^~Pn i+1 with R S {ni, m + 1, . . . , n i+ i\. 
ti H 
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In Remark 4 and in what follows we omit fixed parameters of functions P, /i, and P' . 

In order to complete the proof of Theorem 1 in the case r = 1, p = 0, it remains to note that by definition 

pl,0 _ spB 2n p/ • 

r R 2 +2R ~ L-m=\ R 2 +2R r m l - e -l 



P R 2 +2R — 7 , D2 , D^n ( 14 ) 



i? 2 + 2R R2+2R ^ R 2 + R 



n=l 



According to Lemma 1 and Remark 3, the limit of the right-hand side of equality (fl"4]) for R — > oo coincides 
with the corresponding value of /i, and the limit of the left-hand side (due to Remark 4) is lim Ph° . 

R— >oo 

One can prove Theorem 1 analogously in the case r = p = 1. According to 3J, one can reduce the general 
proof with r = 1 to these two cases. Unfortunately, partitions of the integration interval [0, 1) onto equal 
parts do not enable us to prove Theorem 1 with r > 1. More intricate techniques are necessary. 

Proof of Theorem 1 in the general case. Let n E N. In order to calculate integral fl}, let us divide 
the segment [0, 1) onto p + r + 2nr parts with left ends at the points x + (q) — n; here q — i — 1 + n(p + nr), i 
stands for the number of an interval. Let us calculate the value of the integrand on the interval at the same 
point x+(q) — n (further x(q) = x+{q)). We obtain 

(n+l)0+(n+l)r)-l 

M= lim V (x(q + l)-x(q)) I(x(q)). (15) 

n — >oo * — » 

q—n(p+nr) 

In accordance with Remark 3 and due to the boundedness of the sums in the right-hand side of (|15p , we 
have 

n(p+nr) N-l , . , . 

„ x( g )-x( g -l) ) x(g ) 

' «^oc. ^ n v y " n^oo ^ x(N) y y " 

q=l q=0 L V /J 

Evidently, in the denominator of the latter expression we can write x(N) in place of Lx(A^)J . As a result, we 
get the limit in form (fl~3|) . where tj = x(i) — x(i — 1). It is well known that (see [26l Chap. 3, theorem 14]) 

t\I\ + . . . + tpflN a/ t^Il + ■ ■ ■ + %/iV 

JN = > 1^ v=^ N — > A*' 

tl + . . . + tjV II + . . . + tjv 

provided that the following conditions hold: 

AT AT 

tg/ijy < const y^t' q /t' N . 

q=l q=l 

Let us verify them in the case i 9 = x(q) — x(q — 1), t' q = 1. The first condition x(g — 1) + .x(g + 1) < 
2a;(g) follows from the convexity of the function x(q), the second one does from the existence of the limit 

lim N (l - ^p^ ) = 1/2. Therefore, we get 



N 



N-l R 

» = J im E '(sfa))/^ = lim J2l(x(q))/R. 

q=0 9=1 

By definition Ylf=i I{ x (l))/R = Pr P ■ Theorem 1 is proved. 

3 The gradual "Nose-Hoover" algorithm (the proof of Theorem 2) 

. . and before he knew what he was doing he lifted up his trunk and hit that fly dead with the end of it. " 

R.Kipling. "The Elephant's Child" 

In what follows, without loss of generality, we consider only continued fractions of positive irrational 
numbers. 



G 



The standard "Nose- Hoover" algorithm for finding a continued fraction of a real number x was proposed 
by Klein in 1895 [25] (the term was introduced by B.N. Delone). It is well known p] that this algorithm 
is reduced to the geometric method which constructs the boundary of the convex shell of the set of integer 
nonnegative points (u, v) located above (below) the straight line v — xu. Let eo stand for the vector (0, 1), 
d = (1,0). Put 

e„+i = e„_i + a„_ie„, n = l,2, (16) 

where a„_i is the maximal integer such that the vector e n+ i lies below (for odd n) or above (for even n) 
the straight line v — xu. Klein noted that et„, n = 0, 1, . . ., coincide with partial quotients of the continued 
fraction of the number x. The geometric algorithm results in two polylines which represent parts of sails of 
the continued fraction. 

We need a slightly modified algorithm which results in one infinite polyline L, originating at the point (1, 1), 

whose segments are the vectors a„_ie„, n = 1, 2, Evidently, e n +i + e„ = e„ + e n +i- Consequently, if in 

the standard "Nose-Hoover" algorithm (fl"6|) we add the extra vector e„ and thus go out of the line v = xu 
(i. e., we consider the vector e n -i + (a n -i + l)e«), then we get the first integer point on the segment of the 
polyline of another sail constructed at step n + 1. This fact justifies a simple geometric algorithm which 
constructs the polyline L. 

We begin the construction process with the point (1, 1); it is convenient to connect it with the origin of 
coordinates by the segment which does not enter in L. As the main "constructive" term at the first (nth) 
step we choose the vector e\ (the vector e„). We add this vector till we go out of the line v — xu. The newly 
added vector e ra +i is directed from the origin of coordinates to the point obtained as a result of the latter 
addition up to the step out of the line. See Fig. 1 for the first segments of the polyline for x = 

For convenience, we assume (if the contrary is not specified) that the zero partial quotient (i.e., \x\) 
differs from zero. In order to provide the correspondence to the indices of partial quotients, it is also 
convenient to begin the numeration of segments of the polyline L with zero. 

Let us first consider the approximation of a number x by mediants. The study of series of mediants 
is usually connected with the names of the German mathematician Stern, the French watchmaker Brocot 
(he described this series and applied it in clock manufacturing) or (more often) the English geologist Farey. 
The latter paid attention to the "curious property of vulgar fractions" (this was the title of a short letter of 
Farey published in "Philosophical Magazine" in 1816). Actually, the Farey sequence was described by Haros 
(see [27l p. 36-37]); he considered its application in 1802. Even ancient Greek mathematicians were able to 
enumerate rational numbers with the help of mediants. 

Let fx = /2 = 2j| be irreducible fractions such that Vi , u\ , v% , u-x > and /i < f 2 . Let us define 
the operation \ of finding the mediant (the "insertion" operation) by the formula f\ \. f% = where 
v = vi + v 2 , u = u\ + U2- Earlier we used the denotation introduced by A. A. Kirillov. Evidently, in the 
geometric representation of the fraction - as the integer vector (u, v) the operation J, corresponds to the 
addition of vectors (ui,v±) and (u 2 , v 2 ). The obtained diagonal of the parallelogram appears to be "inserted" 
between its sides. 

Let x € (/i, /2), = /i I f 2 . Since /3 S (/i, f 2 ), one can treat as an approximation of the number x. 
If x G (/i, /s), then we put fi = fi I fs, otherwise we do f± = /a J, f 2 . The process of the approximation of 
the number x by mediants consists in the repetition of these operations for the corresponding intervals. 

The condition x > means that x £ ( j. The initial approximation of x for such an interval is the 
fraction y. One can easily see that for the mapping 

v 

the fraction — the point (u, v) 

u 

the geometric representation of the algorithm for approximating the number x by mediants is completely 
identical to the algorithm for constructing the polyline L. The technical distinction consists in the following 
fact: earlier we constructed each segment of the polyline L "at once", but now it "grows gradually" due 
to the stepwise addition of the next vector &i. We call a simplified version of this algorithm as applied to 
quadratic irrationalities the gradual "Nose-Hoover" algorithm. 

Let (u n ,v n ) stand for the nth integer point on the polyline L, n = 0, 1,2, . . .. We obtain it at the nth 
step of the gradual "Nose-Hoover" algorithm. In Fig. 1, for example, {uq.vo) = (1,1), (ui,i>i) = (1,2), 
(u 2 ,v 2 ) = (2,3), (w3,i>3) = (3,4), (u4,v 4 ) — (5,7), etc. The number i of the vector which we add at the 
nth step, generally speaking, does not exceed n; the identity i = n takes place only for the golden ratio. 



7 




8 



Figure 2: The arithmetic progression rule 



The main inconvenience of the standard "Nose-Hoover" algorithm consists in the following fact: the 
polylines very quickly tend to the straight line u = xv and the geometric illustration requires the high 
accuracy ("all noses become longer"). In the case of a quadratic irrationality we can use only conventional 
illustrations which allow us to "visualize" the topology of polylines L (and therefore the continued fraction) 
without its exact representation on the integer lattice. In essence, these illustrations coincide with the 
approach (described in [23] ) to the visualization of values of a quadratic form on the banks of its river (we 
use the terminology of 23 ). Below we show the geometric juxtaposition. 

Let us draw one more vector which originates from the point (u„, v„); denote it by e' n . The vector e' n is 
defined by the condition e' n + e; + e" = 0, where is the vector which originates from the point (u n ,v n ) 
and goes along the polyline L, e" = -(m„, v n ) (we considered this vector earlier). The collection of vectors 
{e^,ei,e"}, where each one is defined accurate to the multiplier (—1), is a superbasis 23J. This means that 
any pair of these vectors generates the whole integer lattice. The transition from the point (u n ,v n ) to that 
(u n +i, i> ra +i) corresponds to the replacement of one superbasis with another one; the latter differs from the 
initial superbasis only in one of three its elements. Thus, for the polyline L in Fig. 1 the initial superbasis 
is {(0, 1), (1,0), (1, 1)}, the 1st one is {(1, 1), (1, 0), (1, 2)}, the 2nd one is {(1, 2), (1, 1), (2, 3)}, the 3rd one is 
{(1,1), (2, 3), (3, 4)}, etc. 

Let us extend vectors e' n up to the rays which originate at the points (ii„,ij„). As a result, the first 
quadrant appears to be divided onto disjoint connected domains (see Fig. 1). Let us associate each vector 
(u ni v n ) with a domain, whose boundary contains all points of polylines which include the vector (u n ,v n ) 
in the corresponding superbasis. In the picture these domains look like narrow "crevices" located at the 
north-cast of the points (u n ,v n ). We associate the domains which border on the coordinate axes with the 
unit vectors of these axes. 

Without loss of generality, in what follows we assume that in this statement of Theorem 2 the coefficient 
r is positive. 

If, in addition, the equation has a unique positive root, then 

thepoint(u n ,v n ) is located above the line v — x + (q)u rv\ + pv n u n — qu^ > 0. (17) 

If both roots are positive, then this equivalence does not necessarily take place for all n. The quadratic form 
is also positive at the points located below the line v — x_ (q)u. However, since the polyline L tends to the 
straight line v = x + (q)u (it is known that \^ — x\ < l/u^i), correlation (jTTJ) is true for all n > uq. Let Iq 
stand for the number of the segment of the polyline which contains the point (u no ,v no ). In the case of one 
positive root, l = 0. 

Let us write the values of the quadratic form rv 2 +pvu — qu 2 at the points (m„, v n ) in the corresponding 
domain. See Fig. 1 for the result obtained in the case r = 1, p = 0, q — 2 (the upper right corner of the figure 
contains the value calculated at the point (5, 7)). From (fT?)) it follows that with all n > no the points of the 
polyline L belong to boundaries of domains with positive and negative values. The collection of superbases 
which correspond to these points is called the river of the quadratic form [S3] . 

Let us adduce the main result of the theory of quadratic forms which enables one to make the calculation 
of the period of a continued fraction much easier. Using the values of the quadratic form at three vectors 
from any superbasis, one can easily restore the continued fraction. To put it more precisely, the following 
correlation [33J is true. Let domains associated with values a, 6, c, d of the quadratic form be located in 
accordance with Fig. 2. Then values d,a + b,c form an arithmetic progression. Moreover, if h is the common 
difference of this progression, then 

h 2 /4-ab = A/4. (18) 

Therefore, using the values (a, b) for two neighboring domains with different signs separated by the 
polyline L, and the value h, one can unambiguously restore all subsequent values of partial quotients of 
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the number of a segment of the polyline 





1 


1 


2 


2 


a is the value located at the north west of L 


1 


2 


1 


1 


1 


b is the value located at the south east of L 


-1 


-1 


-1 


-2 


-1 


h is the common difference of the arithmetic progresion 


2 





-2 





2 



Table 1 : All possible values of (a, b, h) for the polyline represented in Fig. 1 



the continued fraction, moving along the river of the quadratic form (see the figures in Chapter 1 of book 
[23]). Moreover, according to the arithmetic progression rule, one can restore all previous values of partial 
quotients (whose numbers exceed Iq). 

Remark 5. If h = 0, then the subsequent values and the previous ones coincide. In this case the periodic 
part of the continued fraction is palindromic. In particular, this takes place for p = ; i. e., for the root of a 
rational number (see the latter paragraph of this section). 

Remark 6. The uniqueness of the river of a quadratic form (see \2S$ ) implies that indefinitely restoring the 
previous values, we obtain (beginning with certain n' ) the basic (infinite) part of the continued fraction of 
the associated irrationality X—(q). Hence, one can obtain the period of the continued fraction X—(q) from 
the period a m ( r ,p,g)+i) ■ ■ ■ , a m(r,p,q)+T by operations of the inversion and the cyclic shift. In particular, for 
r = 1, when X—(q) is connected with x+(q) by the evident unimodular transform (it is well known that the 
latter preserves the period of a continued fraction) , we observe the palindromicity. 

Let us now immediately prove the theorem. Its assertion directly follows from equality (118|) . Really, as was 
mentioned above, one can perfectly restore partial quotients dj with i > Iq, using the gradual "Nose-Hoover" 
algorithm with n > uq. This can be done not only with the help of the exact geometric representation in the 
coordinate plane, but also by a "caricature" representation of the motion along the river of the quadratic 
form. Assume that at step n\ we get the same parameters of the arithmetic progression as those obtained 
at step Uq. Let stand for the number of the polyline segment which contains the point (u ni ,v ni ). Then a 
part of the continued fraction a; +i, ai Q +2, ■ • • > a h becomes periodic. In addition, the number l\ — Iq is even, 
because the points (u no ,v no ) and (u ni ,v ni ) are located to one side of the polyline L. In the case of an odd 
period the sequence et/ 0+ i,a; 0+ 2, ...,ai x contains at least two repeating subsequences. In accordance with 
the algorithm we have Yli=i a +i a i — n i ~ n o- Therefore, suffice it to estimate the number of all possible 
triplets (a,b,h) which satisfy (JT5J). The obtained bound obeys formulas ©, (fTU)) . The multiplier 2 in the 
right-hand sides of these formulas appears, because we have to take into account various signs of h, and the 
term r(A/4) does when we consider the case h — 0. Theorem 2 is proved. 

Let us consider the case of one positive (and one negative) root, which is equivalent to rq > 0. If 
|_x+(g)J = 0, i.e., ao = 0, then the numeration of segments of the polyline shifts by one. The above 
considerations imply that only three alternatives are possible. Namely, a continued fraction can have no 
preperiod (m(r,p, q) — 0), or m(r,p, q) = 1, or m(r,p, q) = 2 and ao = 0. In the second and third cases the 
latter partial quotient of the preperiod is less than the latter partial period of the continued fraction. 
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